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AN/SQS-23  Digital  EM  Classifier 
(K.  W.  Harvel,  W.  D.  Howard,  and  F.  L.  Weisser) 


(U-J0UCT)  The  two  weeks  of  submarine  target  data  collection  were  completed 
near  the  end  of  the  current  reporting  period.  The  FM  equipment  was 
transported  to  San  Diego  and  installed  aboard  USS  ROGERS  (DD  876) 
with  a minimum  of  difficulty.  We  were  particularly  pleased  to  be 
able  to  use  spare  sonar  and  fire  control  cables  for  all  of  the  signal 
connections  required,  thereby  avoiding  the  task  of  stringing  cable 
through  the  ship.  Main  power  for  the  EM  equipment  was  tapped  from 
Xhe  fire  control  switchboard. 


(U-JS0W7)  * Defense  Research  Laboratory  (DRL)  personnel  assisted  in 

calibrating  the  receiving  portion  of  the  sonar  system.  Repeating 
the  installation  checkout  required  many  calibration  adjustments, 
the  most  troublesome  of  which  was  the  correction  of  tape  path  align- 
ment errors  on  the  PME. 


1 >1  f 4 


(U Standard  acoustic  measurements  were  made  on  the  system 

transmitter  and  receiver.  Audio  scanner  receive  sensitivity  of 
-17  dB  re  1 .pbarv  was  measured,  with  side  lobe  suppression  a 
relatively  poor  13  dB.  The  receive  bandwidth  was  modified  for  the 
testing  project;  the  measured  response  is  shown  on  Fig.  1.  A few 
individual  stave  measurements  were  performed  with  the  results  shown 
on  Figs.  2 and  3.  Figure  2 shows  the  receiving  response  of  stave  1, 
and  Fig.  3 shows  a comparison  of  the  acoustic  signal  response  of 
several  transducer  staves . The  staves  were  connected  to  the 
preamplifiers  during  these  measurements. 

(U-JQUO)  Source  level  was  observed  to  vary  dep^ding~otT^he~transmiTter- 
combination  used.  With  either  the  "A"  or  "B"  transmit  section  the 
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FIGURE  3 
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source  level  was  measured  to  be  136  dB,  while  with  "A  and  B" 
transmitters  used,  the  source  level  decreased  by  2 l/2  dB.  "A  and  B" 
mode  was  used  for  the  data  recording.  It  was  realized  that  the 
transmission  system  needed  to  be  recalibrated,  but  by  the  time  the 
deficiency  was  noted  the  task  could  not  be  accomplished  before 
scheduled  sailing  time,  and  the  loss  of  source  level  was  estimated 
to  be  of  insufficient  importance  to  limit  our  data  collection. 

Using  the  136  dB  source  level  the  beam  patterns  and  sound  field 
of  the  first  four  harmonic  frequency  components  were  measured.  No 
effort  was  made  to  reposition  the  hydrophone  at  the  higher  frequencies, 
so  the  levels  reported  may  be  less  than  the  maximum  value.  Figure  k 
shows  the  KDT  envelopes  at  the  fundamental,  5 kHz,  and  the  second 
harmonic,  10  kHz.  The  15  kHz  component  showed  no  major  lobe 
structure,  while  the  major  lobe  returned  at  20  kHz.  The  levels  are 
summarized  in  the  table  below: 


Frequency 

Amplitude 

Measured  6 dB 

Component 

re  1 pbar 

Pulse  Length 

5 kHz 

136  dB 

25.3  msec 

10  kHz 

118  dB 

16.7  msec 

15  kHz 

98  dB 

— 

20  kHz 

108  dB 

10  msec  (3  lobes) 

At  20  kHz  there  are  three  lobes  of  near-equal  amplitude,  two  lobes 
being  displaced  28  deg  either  side  of  the  main  lobe,  approximately 
corresponding  to  the  expected  main  lobe  ambiguity  from  elements  spaced 
two  wavelengths  apart. 

The  purpose  of  the  current  data  collection  period  is  to  obtain 
the  recordings  necessary  to  compare  the  classification  performance 
of  several  systems.  ASPECT,  sonar,  M and  "fM-PAIR"  data  were 
emphasized.  Thirty- four  of  the  very  useful  £lpha/Sierra  data 
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5 kHz  - SOURCE  LEVEL  - 136  dB 


10  kHz  - SOURCE  LEVEL  - 118  dB 


FIGURE  4 

RDT  TRANSMISSION  COMPONENTS  (U) 
USS  ROGERS  (DD  876)  AUGUST  1968 
SCALE:  20  msec /div 
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sequences  and  ninety  of  the  Bravo/Sierra  runs  were  completed.  This 
is  almost  twice  the  estimated  minimum  amount  of  submarine  data 
required  to  perform  an  operator  classification  comparison. 

Nonsubmarine  records  were  obtained  for  six  Alpha/November  sequences 
and  seven  Bravo/November  sequences — about  15 % of  the  estimated 
minimum  nonsubmarine  data.  "Chatterbox"  echoes  recorded  at  the 
request  of  Bendix  Corporation  were  obtained  for  twelve  submarine  and 
three  nonsubmarine  sequences . 

(C)  Data  collection  will  resume  in  early  October  for  the  recording 
of  nonsubmarine  segments.  An  additional  requirement  was  added  for 
this  second  recording  period:  sufficient  data  were  needed  to  compare 
operator  detection  performance  with  broadband  vs  narrowband  rever- 
beration background.  Alerted  and  nonalerted  minimum  detectable  signal 
levels  will  be  measured  at  the  AN/SQS-25  playback  facility.  An 
M vs  cw  pulse  audio  experiment,  which  was  reported  previously,  indi- 
cated a performance  improvement  of  8 dB  for  fully  alerted  detection 
in  broadband  vs  narrowband  interference,  respectively.  The  test 
was  not  sufficiently  realistic  to  anticipate  that  the  same  gain  can 
be  realized  in  practice,  and  the  current  data  should  allow  a much 
better  measure  of  relative  operator  performance. 

(C)  Design  of  the  circuit  components  necessary  to  perform  an 

experimental  M modification  to  the  AN/SQS-23  has  continued.  The 
original  transmit  oscillator  was  composed  of  three  crystal  oscillators 
operating  at  l8,  20,  and  22  kHz,  two  discrete  component  flip-flops, 
which  divide  the  oscillator  frequency  to  4.5,  5>  and  5.5  kHz,  a 
transmit  enabling  gate,  a low- pass  filter,  and  a transformer  coupled 
output  stage.  The  redesigned  oscillator  retained  all  the  functions 
of  the  original  oscillator  and,  in  addition,  three  selectable 
JM  oscillators  with  center  frequencies  of  4.5,  5,  and  5.5  kHz  on  an 
original  size  printed  circuit  board.  Operator  selectable  M or 
cw  mode  is  accomplished  by  the  addition  of  a momentary  switch  on 
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the  sonar  control  indicator,  while  the  selection  of  a particular 
frequency,  whether  cw  or  FM,  uses  the  existing  frequency  deviation 
switch  (no  rewiring  is  necessary). 

(c)  To  retain  all  the  functions  of  the  original  oscillator  board,  and 
still  have  the  board  space  necessary  for  the  addition  of  the  FM  capa- 
bility, two  TTL  integrated  circuits  in  the  a.4  pin  dual- in-line 
configuration  replaced  the  discrete  flip-flops  and  the  transmit 
enabling  gate.  The  crystal  oscillators,  the  filter,  and  the  output 
stage  were  retained.  The  FM  oscillator  is  a unijunction  transistor 
relaxation  oscillator  whose  capacitor  charging  rate  is  varied  by  the 
ramp  generator  and  the  transmit  frequency  deviation  switch. 

(c)  The  FM  oscillator  and  the  ramp  generator  are  very  similar 

circuits.  The  capacitor  charging  rate  of  the  ramp  generator  is  held 
constant  by  a constant  current  source,  resulting  in  a linear  ramp. 

The  time  constant  of  this  RC  oscillator  is  chosen  to  obtain  a period 
of  50  msec . The  ramp  generator  is  free  running  and  is  not  synchro- 
nized to  the  enabling  gate  on-time.  The  FM  oscillator  capacitor 
charging  rate  is  changed  by  the  linear  ramp,  resulting  in  a linear 
FM  slide  that  varies  l60  Hz  either  side  of  the  selectable  center 
frequency . 

(c)  The  decision  to  use  a rather  simple  unijunction  FM  oscillator 
instead  of  a more  precise  voltage  controlled  oscillator  was  based 
primarily  on  the  amount  of  board  space  available.  A single  uni- 
junction oscillator  is  capable  of  producing  all  three  center 
frequencies  required  with  a minimum  of  space  consuming  circuitry. 

The  major  disadvantage  of  such  a circuit  is  the  difficulty  of 
achieving  frequency  stability  with  varying  temperature.  In  order 
to  obtain  an  acceptable  variation  of  less  than  20  (±10)  Hz  along 
any  portion  of  the  FM  slide,  it  has  been  necessary  to  temperature 
cycle  each  board  to  select  the  proper  temperature  compc  '.sating 
resistor. 
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(C)  It  should  be  pointed  out  that  an  attempt  was  made  to  achieve 
crystal  controlled  PM  slide  center  frequencies  by  mixing,  but 
although  two  different  mixing  circuits  were  designed,  the  entire 
slide  frequency  was  still  c.ependent  upon  the  stability  of  the  least 
stable  (PM)  oscillator. 


B.  Systems  Analysis 
(S.  P.  Pitt) 


( U- petJCT  Work  continued  in  all  areas  mentioned  in  the  previous  progress 
report.  However,  only  studies  where  specific  information  was 
generated  will  be  mentioned  here. 


1 . Floating  Point  Correlator 

(c)  The  data  necessary  to  complete  this  study  were  reprocessed 

during  this  quarter  after  it  was  discovered  that  the  quantization 

routine  in  the  quadrature  correlation  version  of  the  program  was  not 

precisely  correct.  The  differences  lay  in  the  lower  order  (least 

significant)  bits.  Comparison  of  these  data  after  reprocessing 

showed  only  slight  differences  in  nearly  all  cases . The  correlator 

study  is  now  essentially  complete  and  will  be  described  in  a technical 

memorandum  shortly.  Recently,  however,  the  question  of  how  useful 

this  type  of  quantization  is  for  beamforming  has  been  raised.  The 

possibility  of  quadrature  beamforming  digitally  using  floating  point 

quantization  has  been  discussed  and  partially  analyzed,  the  result 

is  as  follows.  Using  four  bits  (exponential)  quantization,  a dynamic 
15 

range  of  20  log  2 = 90  dB  can  be  achieved.  Using  quadrature 

sampling,  beamforming  can  be  done  at  low  frequencies  in  times  that 
would  apparently  allow  digital  formation  of  many  beams  in  the  time 
between  samples  for  each  beam,  for  typical  sonar  bandwidths,  allowing 
electronic  beam  scanning.  Whether  the  errors  involved  will  allow 
accurate  enough  beamforming  remains  to  be  seen;  a preliminary 
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to  2/w,  where  W is  the  bandwidth  of  the  impinging  signal, 

X(t  - -r)  = X(t),  Y(t  - t)  = Y(t),  so  that  the  transformation  becomes 
a simple  rotation. 

A beam  formed  by  this  type  combination  of  several  receivers 
would  have  the  same  characteristics  as  the  present  beams  formed  in 
the  AN/SQS-23C  preformed  beam  receiver,  except  for  the  errors  due  to 
quantizing.  A short  study  to  determine  the  efficiency  of  this  type 
of  beamforming  has  been  started:  the  complicated  quantizing  schemes 
already  being  used  in  the  correlator  study  should  allow  rapid  simu- 
lation and  testing  of  the  techniques.  A technical  memorandum  is 
anticipated  upon  completion  of  the  study. 

2.  Covariance  Analysis 


The  program  to  produce  covariance  analyses  for  input  signal 
ensembles  has  been  completed  and  used  extensively  in  analysis  of 
ASPECT  data.  The  beam  data  for  the  submarine  reported  previously 
were  so  analyzed  and  plotted  on  a graph  with  abscissa,  the  time 
between  members  of  the  ensemble,  and  ordinate,  the  normalized  cross- 
correlation coefficient,  as  shown  in  Fig.  5.  The  periodic  nature  of 
these  data  is  quite  obvious.  This  analysis  led  to  the  digitizing  and 
analysis  of  echoes  from  other  aspect  angles.  The  analysis  is  nowhere 
near  complete,  but  preliminary  results  indicate  that  the  same  perio- 
dicity exists  for  stern  aspect  data.  The  most  reasonable  explanation 
for  the  periodicity  is  the  oscillation  of  the  submarine  about  its 
desired  course.  Naval  Ship  Research  Development  Center  was  contacted 
to  obtain  data  which  would  support  this  as  a possibility;  expected 
amplitude  of  such  oscillations  and  data  have  been  promised.  In  order 
to  demonstrate  that  this  is  the  true  explanation,  it  is  planned  that 
data  will  be  analyzed  for  many  different  conditions  on  the  same 
target,  i.e.,  depth,  speed,  and  aspect  angle.  By  correlating  the 
results  with  expected  results  using  the  "oscillation"  hypothesis, 
its  validity  (hopefully)  can  be  determined. 
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3.  Multipath 

(U-B0WT)  The  problem  of  estimating  and  eliminating  the  effects  of 

multipath  from  sonar  echoes  has  been  pursued  following  examination 
of  these  data,  described  in  the  previous  progress  report.  An 
analytical  treatment  of  the  problem  follows: 

Consider  a signal 


where 


X(t)  = f(t)  + a f(t  - T)  , (1) 

a 

f(t)  is  the  signal  to  be  expected  in  the  absense  of  multipath,  and 
T is  the  delay  of  the  second  path. 


Let  a ± 1;  for  convenience,  a < 1.  There  are  several  ways  to 
represent  the  recovery  of  f(t)  from  X(t),  T and  a being  known  or 
estimated. 

In  the  frequency  domain, 


X(f)  = f(f){l  + a e"2jtlf]) 


from  which 


f(f)  = 


1 + a e 


with  f(t)  being  recovered  according  to 


f(t)  s J f(f)  e2nlft  df 
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(U-S9W}  With  X(  f ) being  analytic  in  the  upper  half-plane, 

we  might  close  the  contour  above  and  evaluate  the  residues  of 


lill 


2itift 


1 + a e 


■2jtifT 


at  the  zeros,  summing  then  to  get  f(t). 


(U- MdUU)'”-  However,  there  are  far  more  convenient  inversions.  Since 
a < 0,  we  expand  Eq.  (3)  for  real  f: 


A A - 

F(f)  = X(f){l  - a e 
and  invert  to  get 

f(t)  = X(t)  - aX(t  - T)  + a2X(t  - 2T)  - a3X(t  - 3T)  . . . 


-2jtifT  2 -2jtif2T  3 -2*if3T 


+ a e 


a e 


• • •}  > (5) 


= x(t)  + Y (-a)n  X(t  - nT) 


(6) 


This  agrees  with  intuitive  notions  for  recovery  of  f(t)  from  X(t). 


( U- B0U0 Equation  (6)  may  be  obtained  also  by  using  an  operator 

notation,  introducing  a time- translation  operator  S(t)  = eTd/ot^ 


Then 


f(t  - T)  = 

S(-T)f(t)  . 

(7) 

X(t)  = 

[1  + aS(-T)Jf(t) 

(8) 

f(t)  = 

[1  + aS(-T)]"1  X(t ) ; 

(9) 
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The  operator  maybe  expanded: 


[1  + aS(-T)]'1  = 1 - aS(-T)  + aS(-T)  X aS(-T)  + . . 


= 1 -aS(-T)  + a2S(-2T)  - a5S(-3T)  + . 


00 

f(t)  = X(t ) + ^ (-a)n  S(-nT)  X(t) 

- l" 


which  is  equivalent  to  Eq.  (6). 

(U-P0WC7  A simpler  implementation  of  the  elimination  process  can  be 

obtained  for  signals  of  finite  length  by  assuming  that  epoch  (t ) 
can  be  determined  and  that  the  leading  edge  (T  sec)  of  the  signal  X(t) 
is  purely  f(t).  Then 


X(t  + t)  - f(t  + t) 

f(t  + t - T)  = 2 2 = 0 , 0 = t < T , (12) 


so  that 


f(tQ  + t)  = X(tQ  + t)  - af(tQ  + t - T)  , T = t < 2T  . (13) 


The  same  process  is  repeated  for  2T  = t < 3T,  using  the  reconstructed 
f(tQ  + t)  in  the  last  term.  Hence,  by  storing  the  entire  signal,  one 
can  obtain  the  delayed  version  of  the  signal  by  a simple  subtraction 
process.  The  problem  remaining  is,  of  course,  that  a and  T must  be 
known  or  estimated. 


(U-«JUUT  A computer  program  to  implement  Eq.  (13)  was  generated  during 

this  quarter.  The  values  of  a and  T were  estimated  in  two  ways, 
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(l)  by  visual  inspection,  and  (2)  by  use  of  the  autocorrelation 
function  for  X(t).  It  is  obvious  by  inspection  that,  given  X(t)  as 
in  Eq.  (l),  the  autocorrelation  function  t)  = E[X(t)  X(t  - t)] 
will  contain  relative  peaks  at  (at  least)  % = 0 and  t = T.  A 
computer  program  was  written  to  compute  0(  t)  and,  from  0(t),  estimate 
a and  T.  The  program  has  been  tested  on  analytically  generated  data 
and  real  data.  The  results  of  the  study  will  be  reported  in  the  next 
quarter. 

4.  Digitized  Sonar  Stave  Data 

During  this  quarter  Professor  E.  A.  Patrick  and  Mr.  L.  Shen 
visited  DKL  for  a cooperative  program  between  DRL  and  Purdue 
University.  Both  DRL  and  Purdue  initiated  programs  to  examine  sonar 
data  extracted  by  both  space  and  time  sampling  of  the  acoustic  field. 
At  present,  DRL's  effort  is  directed  primarily  to  an  examination  of 
the  reverberation  and  noise  fields  and  beamforming,  whereas  Purdue  is 
studying  the  representation  of  echoes  for  recognition  purposes.  For 
both  groups,  sampled  data  from  the  output  of  individual  staves  of  the 
AN/SQS-23  sonar  system  were  required.  The  programming  necessary  to 
obtain  such  samples  was  generated  at  DRL  and  these  data  were  obtained 
as  follows . The  Ampex  FR-l800  was  used  to  record  the  outputs  of  the 
ten  staves  in  the  direction  of  the  target  during  playback  of  sonar 
tapes  on  the  SME.  These  data  were  recorded  at  60  ips  along  with  the 
reference  track  from  the  SME  tapes  and  the  output  of  the  audio 
beamformer.  By  careful  manipulation,  all  ten  channels  were  digitized 
in  quadrature  for  most  of  these  data  recorded,  which  consisted  of 
submarine  echoes  at  various  aspect  angles.  Copies  of  these  data 
were  forwarded  to  Purdue  University.  Time  has  not  permitted  a 
significant  amount  of  analysis  by  either  group.  Nonsubmarine  data 
must  eventually  be  obtained  in  the  same  manner  for  further  development , 
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(U-FOUCr)-  One  interesting  observation  was  made  during  handling  of 

these  data;  by  comparing  the  signal-to-noise  ratio  (s/N)  on  indi- 
vidual staves  with  that  from  the  audio  beamformer,  it  was  found 
that,  at  least  for  one  case,  the  gain  achieved  by  the  beamformer  was 
about  12  to  13  dB.  The  theoretical  gain,  assuming  independent 
Gaussian  noise  at  the  staves,  is  about  1 6 dB.  The  difference  could 
be  caused  by  dependence  of  the  noise  samples,  by  incoherence  of  the 
echo  from  stave  to  stave,  by  inaccurate  beamforming,  or  by  all  three. 
All  of  the  possible  causes  are  being  examined. 

5 . Clipped  Correlator 

( U- P8UQ-)~  The  analysis  of  the  inconsistent  output  of  the  clipped 

correlator,  discussed  in  the  previous  progress  report,  was  post- 
poned because  of  lack  of  availability  of  the  A/D  and  because  of 
the  sea  trip.  It  is  expected  that  the  study  will  be  completed 
during  the  next  quarter. 

6.  Comparison  Between  Direct  Quadrature  Correlation  and 
Quadrature  Correlation  Using  FIT 

( U- P0UO  Y'  Experiments  with  real  data  were  run  during  this  quarter  to 

compare  correlations  using  the  direct  method  with  those  using  FIT. 

A program  that  computes  the  Fast  Quadrature  Transform  was  developed 
and  used  by  the  Computer  Science  Division  to  generate  envelopes  of 
correlation  functions  from  quadrature  sampled  data.  The  machine 
language,  direct  quadrature  correlator  program  of  Systems  Analysis 
was  used  on  these  same  data.  The  outputs  are  compared  in  Fig.  6, 
where  it  is  seen  that  they  are  identical.  For  this  case,  where  a 
maximum  of  ^12  data  points  were  used,  the  direct  method  was  faster. 

It  is  expected  that  for  more  points  the  FQIT  method  could  probably 
be  streamlined  (using  machine  language)  to  make  it  somewhat  faster, 
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whereas  the  direct  method  has  been  optimized  in  speed  as  much  as 
possible  on  the  CDC  3200. 

( U itjHJO )'  The  paper  on  quadrature  sampling  will,  according  to  the 

present  schedule,  appear  in  the  correspondence  section  of  the 
November  issue  of  The  Journal  of  the  Acoustical  Society  of  America, 
instead  of  the  October  issue,  as  reported  earlier. 
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